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13. ABSTRACT 

This program uses Lemke, Salkin and Spielberg's (ref. 1) Set Covering Algorithm 
(SCA) to optimize a traffic model problem in the Scheduling Algorithm for Mission 
Planning and Logistics Evaluation (SAMPLE). SCA forms a submodule of SAMPLE and 
provides for input and output, subroutines, and an interactive feature for per- 
forming the optimization and arranging the results in a readily understandable 
form for output. 
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SET COVERING ALGORITHM, • 

A SUBPROGRAM OF THE SCHEDULING ALGORITHM FOR 
MISSION PLANNING AND LOGISTICS EVALUATION 

- • r 

INTRODUCTION 

This documentation provides a general deseription of the Set Covering 
Algorithm (SCA) computer program, and includes functional specifications, 
functional design and flow, and a discussion of the program logic. The SCA 
program is a submodule of the Scheduling Algorithm for Mission Planning and 
Logistics Evaluation program (SAMPLE) and has been designed as a continuation 
of the Mission Payloads (MPLS). The MPLS uses input payload data to form a 
set of feasible combinations which are collections of payloads that meet cer- 
tain system constraints (e.g.. Shuttle weight- to-orbit capability); from this 
combination set the SCA selects a subset with minimum cost such that all pay- 
loads are contained without redundancy. The subset of feasible combinations 
is called a traffic model. To date, the program has had two main uses: 

a. To provide input data for the Operations Simulation and Resource Scheduling 
(OSARS) submodule, and 

b. To provide the user a tutorial option so that he can choose an alternate 
traffic model in case a particular traffic model cannot be scheduled by 
the OSARS. 

The SCA program was begun in 1974 with an input of less than 50 payloads per 
year and the assumption of unity cost coefficients in the objective function. 
Presently, this program has been expanded to solve 100 payloads per year with 
an option of five different performance criteria in the objective function. 
Since a general SCA program solves a problem with Inequality constraints, an 
appropriate change has been made to convert it to solve the traffic model 
problem, which has aTT equality constraints. 



2.1 DEFINITIONS AND 

2.1.1 DEFINITIONS 

Basis 

Canceled variable 
Ceiling 

Decision variable 
Dominated column 

Extreme point 
Fractional variable 

Free variable 
Level 

Minimum-cost-per'- 

constraint-satisfied 

rule 

Occurrence table 
Partial solution 


2. PROGRAM DESCRIPTION 
SYMBOLS 


- A basis for a n-dimensional euclidean space, 

is a linearly independent subset of vectors from 
E*^ which spans the entire space 

- A variable is said to be canceled if its value is 
set to be zero 

^ The lowest upper bound 

- A zero-one variable which corresponds to a feasible 
mission 

- In a matrix A or TA^-} d - 1, n A^. is said to 
be dominated by A^. if A^. £ A^. for i f j. A. 
is called a dominated column of A. and A. is 
the dominating column of A.. 

J 

- The corner points of a convex polyhedron. 

- The nonintegral variable in the linear programming 
solution 

- A variable which has not been fixed at any value 

- The number of variables fixed at one 

- In the sub-subprobTem, for each variable the value 
of its cost coefficient divided by the number of 
constraints is calculated. The variable with the 
minimum of such value is set to one first. 

- A table for each payload of all feasible combi- 
nations that include that payload 

- A set of variables fixed at one or zero 
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Preferred row 
Preferred set 
Slack variable 


Subproblem 


Sub~subproblem 


Upper bound 


- The row containing the smallest number cf ones in 
the current subproblem's constraint matrix 

- The set of columns in the constraint matrix which 
contains the entry of one in the preferred row 

- In general, it is desirable to convert any in- 
equalities in the constraints into equations which 
are much more convenient to work with in the linear 
programming problem. The conversion is carried out 
by introducing some additional variables which are 
called slack variables. 

- The problem contains only free variables with 
certain constraints deleted from the main problem 
by the satisfaction of the partial solution 

- The problem contains the columns of A associated 
with the positive fractional variables and the rows 
corresponding to the cons'f'raints they satisfy 

- A known value that the solution of the Set Covering 
Problem will not exceed 


2 . 1.2 

A 

c 

e 

F 

I 

i 

a 


SYMBOLS 

An mxn matrix with zeros and ones as elements; constraint matrix 

An mxl matrix of all positive numbers; cost coefficient vector 

An mxl matrix of all ones 

Indices of free variables 

An mxm identity matrix 

The row index of matrix A 

The column index of matrix A 

Level counter 
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M A »!cnnsgative number in the objective function, set to be 

very large to solve a linear programming problem with equality 
constraints; in this program, set M 500 

m Number of constraints in the set covering problem 

n Number of decision variables in the set covering problem 

s An mx1 matrix of slack variables 


y The column vector of zero-one variables 

z Current value of objective functional, i.e,, c‘y 

z* Upper bound for c’y 

a = 1(0) Reaching a point on a forward (backward) step 

Y Iteration counter 

n(il,j) The number of constraints of the level t subproblem 

variable j can satisfy 

^(X,) Lower bound for the minimal functional of level I subproblem 

a Sequence vector, recording nonzero components of y 

T Tolerance parameter about z^ 

(j)(Jl) Number of constraints in the current (level a) subproblem 



2.2 GENERAL DESCRIPTION 


2.2.1 PROGRAM CAPABILITIES 

The SCA program optimizes a traffic model problem over an objective function 
which consists of a set of user-selected performance criteria such as Orbital 
Maneuvering System (OMS) weight, load factor, or Shuttle cargo bay utilization 
for each feasible combination. This problem also includes a set of constraints 
which assures no redundancy of payloads in the traffic model. For application, 
a set covering algorithm developed by C. E. Lemke, H, M. Salkin, and 
K. Spielberg (ref. 1) has been adopted. The main advantage of this algorithm 
is that it permits a rather efficient and simple solution procedure that is 
basically a (zero, one) branch and bound search logic coupled with linear 
programming (LP) and suboptimization techniques. The suboptimization tech- 
nique can construct very good integer solutions from the solutions to LP 
subproblems. 

The formulation of a traffic model problem can be exactly fitted into the 
mathematical model of a Set Covering Problem (SCP), Before a formulation 
example is presented, a general understanding of the form of this model will 
be helpful. 

The SCA solves the SCP which has the form: 

min 

subject to 
and 

When this model is applied to the traffic model problem, is taken as a 
decision variable on a particular feasible combination j (or flight j). 
Flight j is considered to be chosen when x^ = 1, otherwise x. - 0. For 
each payload, there is a correspondent constraint which insures ihe nonredun- 
dancy of this payload in the traffic model . Vector c stores the performance 
criteria for all of the feasible missions; M is set to a large value to 
assure that constraints are satisfied. The application of SCA to the traffic 
model can be easily demonstrated by the following example. , 


c'x + Me's 

Ax - Is = e, s ^ 0 

X, = 0 or 1 
J 


Assume that three payloads have to be launched in a particular month. After 
checking all of the possible combinations, only five are considered to be 
candidates. These candidate flights are called feasible combinations; a 
summary of the payloads they carry is given in the following table. 


Feasible combination (j) Payloads 

1 No, A 

2 No. A, No. B 

3 No. B, No. C 

4 No. B 

5 No. C 


Assume that associated with each feasible combination j, there is a 
cost factor c^. The problem is to formulate a mathematical model for 
determining the traffic model which gives the minimum cost. In this 
example, c. is set to unity, which implies that we are seeking a traffic 
model which consists of the minimum number of feasible combinations. Let 
X. (j = 1 , 2, 3, 4, 5) be the decision variable over the selection of a 

J . 

particular feasible combination j in the traffic model, A unity value 
of X. implies that feasible combination j has been selected to be a 

J 

member of the traffic model. Since cost has been chosen as a measure of 
effectiveness, the object is to minimize 

Z = X, + Xg + X3 + + Xg 

subject to the restrictions developed in the following paragraphs. 

The constraint in this situation is that the same payload cannot be contairied 
in more than one feasible combination in the traffic model. The mathe- 
matical statements of the restrictions for three payloads are 

A: *1 ’‘2 ° 
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B : ^2 *^.^3 ^ ^4 ~ ^ 


Finally, there are the binary restrictions, i.e., x. = 0 or 1 . Therefore, 

iD 

in sunmary, the mathematical model for this problem is the following. Minimize 

z = + X2 + X3 + + X5 

subject to: 




= 1 


*2 + X3 + X4 =1 


X3 + Xj = 1 


and X . = 0 or 1 5 j = 1 , 2, 3, 4, 5 

J 

This formulation can also be written in a matrix form: 


min c'x + Me's 

subject, to Ax - Is = e, s > 0 


X. = 0 or 1 

a 


where 


C = [1, 1, 1, 1, 1] x' = Xg* X3, X^, Xg] 
e‘ = [l, 1, 1] s' = [s^, Sg, S3] 
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1 1 0 0 0 


1 

o 

o 

A - 

0 1. 1 1 0 

I. = 

o 

o 


0 0 10 1 


0 0 1 




L J 


The three feasible solutions are x' ="[1, 0, 1, 0, 0], [0, 1, 0, 0, 1], and 
[1 « 0, Oi T], but the optimum Is either of the first two as they give the 
minimum number of feasible Gombinations. In this simple example* feasible 
solutions can be easily noted by observation, but in dealing with a large- . 
size problem of 100 payloads and 500 combinations, there are 2^^^ possible 
solutions, so a more efficient approach such as the SCA must be employed. 

The details of this algorithm are given in the technical description. 

2.2.2 OPERATIONAL CAPABILITIES 

The operational capabilities of the SCA program were designed to permit the 
user to specify his particular optimization problem according to his needs. 
One of the features is the user's selection of performance criteria over ob- 
jective function. For example, if the user Tikes to see a traffic model of a 
minimum number of missions, he will have the choice of using unity as a per- 
formance criterion, or he may use the surplus Shuttle bay length as a perform- 
ance criterion to find a traffic model which gives the maximum utilization of 
space. At present, the SCA permits the user selection of one of the 
following performance criteria: 

a. Unity minimum number of missions 

b. One minus load factpr; maximum utilization of Shuttle cargo weight 

allowance' 

c. QMS weight: the minimum OMS weight 

d. Payload margin; minirnum payload margin 

e. One minus payload length: maximum utilization of Shuttle cargo bay 

This option can be even more selective if the user has other criteria to be 
defined in the future. Another available option is the suppression or speci- 
fication of certain feasible combinations in the traffic model; the program 
gives the user an error message if he attempts to suppress or specify combi- 
nations which would lead to an infeasible solution that is not a traffic 
model . 
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2.3 TECHNICAL DESCRIPTION 
2.3.1 ANALYSIS 

The analysis for this program is quite lengthy. For a thorough discussion of 
the Set Covering Algorithm, the reader is referred to C. E, Lemke, 

H, M, Sal kin, and K. Spielberg's "Set Covering by Single-Branch Enumeration 
v/ith Linear Programming Subproblems," Operations Research VoT. 19, pp. 998- 
1022, 1971 (ref. 1). 


As this model is formulated, the set covering problem is to minimize c’y 
constrained by Ay > e, yie{0,T>; je{l,2,..n}, which will be referred as (I), 

A = is an m by n matrix with a|j - 0 or 1 , arid e is an m by 1 

column of all I's. Here c is the cost coefficients which are to be positive 

numbers.' 

In the application to the traffic model problem, a small variation about the 
. formulation should be realized. Assuming that 

(a) A row of A corresponds to a payload 

(b) A column of A corresponds to a possible assionment for a feasible 

combination 

(c) The ones in the column denote payloads that could be handled by such an 
assignment 

(d) in a particular solution y, variable yi: = 1 corresponds to actually 

■ J 

using such a assignment, whereas y. = 0 means that such an assignment 

J 

: will /not; be -used,. :■ 

one will note that the constraints Ay e do not require the solution of 
(I) to be a trafflG model , so y^e are really interested in solving a set 
covering problem with equal ity cohstraints, i.e. , Ay - which is usually 
/ called a set partition problem. However, with a Tittle modification in the ' 
formulation, a set partition problem can be solved by SCA without any change 
of the algorithm's generality. 
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Before we discuss the algorithm, certain facts about the development of SCA 
and the application of SGA to the set partition problem will be introduced, 

y*s 

Compared to the problem (I), an LP problem is defined as 


(I^p): min (c'x: Ax > e, x ^ 0) 


1 . 


2 . 


3- 

4. 


5, 


(I) Is feasible if and only if is feasible. (Here feasible 

means that a feasible solution exists). 


Assume that Z 9 snd Zj^p are the optimal solution values of (I) and 


(Ij^p), respectively; then Zj^p is always a lower bound of Zj» V'®*» 


7*'^ > Z* 
i- I ^ ^|_P 


If (I|_p) is solved and integral feasible, then (I) is solved and 


2* = Z 


* 

LP' 


If (I| n) does not give an optimal integer solution, i.e., some x. of 

x^p are not integral, then a rounded-up solution can always be obtained 

by setting all the nonintegral variables to 1 (y. = 1 whenever x. > 1). 

j J 

This rounded-up solution is {I)-feasible and an upper-bound of (1). 


A rounded-up solution obtained from a non-integral extreme point x of 
(I|^p) can always be reduced to another (I)-feasible solution with a 
smaller cost. We have Ax> e and 0 < x ^ e. Consider the columns of 
A associated with the positive nonintegral variables, and the rows 
of A corresponding to the constraints that these variables explicitly 
satisfy. Identify the part of fractional x and the corresponding 
matrix selected as above by a superscript *. Then we have A*x* ^ e and 
0 < X* e; i.e. , every row of A* has at least two 1 's. Therefore, 
setting some of the x* variables to 1 and maintaining A*x* > e gives 
a better integer solution than the rounded-up solution which calls fqr 
setting all x* variables to 1. This reduced integer solution is called 
a purified solution which is reached by successively setting fractional 
variables to I by the '*niinimum-cQSt-per-eoh3tra1nt-sat1sf1ed" rule until 
all variables are either 0 dr 1 and all constraints are satisfied. 



The procedure tp reduce the roundsd-up solution to a purified solution is 
called "purification." 

A' ■ ' ^ ■ 

6/ Let y be any feasible splut ion to (I). If y is not an extreme point 
for (I, i.e. , the column of A corresponding to y, = 1 and the columns 
of -I corresponding to s-> 1 from a linear dependent set, it can be 
reduced to a feasible solution y* for (I), which is an extreme point 

A 

for (I^p) and yields a better value for Zj, We demonstrate this because 
the puri f i ed sol ut i on may not be an extreme point . Let y be any (I ) - 
feasible solution that is not an extreme point for Since 

Ay - Is = e, some of the slack variables must be positive; otherwise, 

Ay - e, which implies y is an extreme point for {Ij^p} since the columns 
of A corresponding to yj " f will be linearly Independent. Permuting 
rows to get positive slack variables last, one obtains 




with s = jA./s^ > l|. Then a permutation of columns to get positive y^'s 
first, leads to 



The last expression means that (1 ) A 2 ^^ - s - e (with s ^ e), i.e. , each 
row of Ag,i contains at least two I’s and U) A-j^e = e, i.e. , each row of 
A^^ contains exactly one 1. Thus, permuting rows and columns further, one 
may exhibit 





with the above, relation (i) may be rewritten as: 



Now the columns of 



clearly form a linear independent set, whereas those of 



2-n 



are linearly dependent by the hypothesis that y is not an extreine point. 

p 2 2 

Therefore, A07 has at least one column, say column A- Then s >e> A. 

Ci p J ^ - *^2 

demonstrates that deleting the A^ column and replacing s by s' = s - A.> 
yields a new feasible solution, with cost reduced by the cost of the deleted 
column. 

Now suppose the above procedure is repeated. Then one is either again not at 
an extreme point and may obtain another cost reduction, or one is at an ex- 
treme point by virtue Of either (a) s = 0, or (b) A^i - 0. 

As an example, consider c* = (1 ,2,1 ,1 ,2,3,1 ) and 


(column) 1 



1 

0 

1 

1 

G 


2 3 4 5 6 7’ 
0 1 0 0 0 1 
110 0 1 0 
0 0 0 0 0 1 
0 0 1 1 10 
10 1 0 0 0 


y {1 ,1 ,0,1 ,1 ,0,0) is (I) feasible with c'y = 6 and s = (0,0, 0,2,1). 
Hence 


(co iumn) 

1 

2 

4 

5 

.3 

6 

7 

r 

A^l 


“i 

0 

0 

1 

0 

0 

0 

0 

|1 
i T 

0 

1 

1' 

0 

^^21 

^22/ 

L 


-L 

_0_ 


0 

1 



1 

0 

1 

1 

■O' 

1 

0 



_o 

1 

1 

0 

1 0 

0 

0 
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*<* 




• 

1 

o 

o 

0 

0 

1 {0 

0 

1 

1 



1 

o 

o 

0 


(coluirm) 12 4 5 





MM . 

1 

0 ! 0 

1 

o" 


0 ■ 


0 

1 ! 0 
» 

0 

a’’ 

A^l 

0 


L 

1 o> 

1 — f 

1 o* 

1 ’ 1 

J1 

a"* 

^21 

A^ 

^21 


1 

0 

0 |i 

1 { 1 

1 

0 

PM. 


Thereforej set = Oj then ” (^) ~ (0" ~ 0,1,0) 

y'= (1 ,1 ,0,0,1 ,0,0) and c’y’= 5. Repeating the procedure, we have 


column 12 5 


A. 


11 


0 

0 


21 ^21 


1 1 
2 0 

3 1 
5 0 

4 1 


0 { 0 

1 i 0 

0 i 0 

1 lo 

-oir 


Finally, set y^ = 0v then s” = 1-1*0 or s" = (0,0, 0,0^0), and 
’V = (1 ,1 ,0,0,0) is an extreme point for (Ij^p) and c‘y" = 3 > c'y = 6. 
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7. Noticing the payload nonrecundancy requirement in the traffic model, we 
now consider a set partition problem (Ay = e) in its relation to the cur- 

yS 

rent problem (I). The actual slack cost is, of course, zero. However, 
by assigning a "high penalty" cost M into all the slack variables, one 
will tend to get a minimal (I) solution with very few positive slack 
variables. For M large enough, we can either find the minimal solution 
or show that it is infeasible. Thus, the equality problem is equivalent 
to 


(I)^: min (c'y + Me's; Ay - Is = e, y. = 0, 1; M large) 

J J 

or a zero-slack cost problem by multiplying the rows of Ay - Is - e == 0 

£ 

by M and adding the sum to c'y, or, (I) is equivalent to 

(I)^ min (c‘y): Ay - Is = e, y. = 0,1} - Me'e (M large) 

y j 

Here c = c' + Me*A; i.e, , to each cost element c. one adds M times 

J 

the number of I's in the associated A columns. 

As an example, consider 

(I.|) iTjir, + 4yg + ly^ + 2y^ 

subject to y^ 1 ^ > 

y^ + V3 1 U 

Vl +V3 

y^, Y 2 > y 3 > y^ = 0 or 1 

Then the minimal y is (0,0, 1,1) and c'y = 3 with s = (0,0,1), Consider 
the problem with equality constraints, and set M = c'e - 12, (I.j) becomes 
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(I^) ; min 5y^ + Ayg + ly^ 2y^ + 12s^ + ^^^2 ^ ^^^3 


subject to y 


1 


yg + y3 

^1 ’ ^ 2 ’ ^ 3 ’ ^4 


" s. 


- s. 


- s. 


0 or 1 ; s^ , Sg, s^ 


= 1 
= 1 
= 1 

= 0 or 1 


After each row of constraints is multiplied by 12 and added to the objective 
function, we get 


min 17y^ + 16y2 + + 26y^ - 36 

subject to y, ■** 

’ ^2 ^ ^3 , -=2 =1 

h * ^3 ^ h '®3 " 

y 2 » y 3 > y 4 = 0 or 1; s^, s^, = 0 or 1 

y\ 

The minimal y is (0,1, 0,1) and c'y 36 = 6 with s=(0,0,0). The above 
procedure permits the solution of the Set Partition Problem (SPP) by means of 
an algorithm (SCA) oriented toward the set covering problem. It is recom- 
mended that M be set to 500 in the traffic model problem. 

2.3.2 METHOD OF SOLUTION 

/% /N 

To solve problem (I), namely, (I); min (c'y: Ay ^ e, y. £ {0,1}, 

* o 

j e {1 ,2,. . . ,n}), an enumerative single-branch scheme is employed. 
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The search starts at the origin (node 0), with all y. ''free", i.e., tenta- 
tively considered to be either 0 or 1 and possibly to be fixed at value 1 on 
a forward step. At the general node v (on "level" fi,), a components of y 
have been explicitly fixed at 1 on forward steps. Others may have been 
"canceled" (fixed at 0) at level A, as one has ascertained that correspond- 
ing forward steps would not lead to a solution. 

The task at node v is to: 

(i) Solve a basic problem U^p) with fixed variables substituted, and to 
look for improved solutions as outlined in 2.3.1. If the objective 
function of (I|_p) exceeds an available bound on (I), or (I^p) is 
infeasible, one may "backtrack", i.e., the search reverts to the 
predecessor node v", linked to v by the branch j* (variable 
y^* having been fixed at 1) and the search continues at node v" 
with y^* canceled to 0. (If no predecessor exists, the search 
terminates and the solution is the existing upper bound.) 

(ii) Cancel whatever variables may be canceled from further consideration. 

(Hi) Select among the remaining variables a "branch" variable j* to be 
fixed at 1 on the next forward step. The state of search is 
essentially recorded in an + 1, n) n matrix: ri(5-,j) = 0 when 

variable j was fixed (i.e., canceled or selected as a branch) at 
some level up to and including the current level H (i.e., at some 
predecessor node of v, or at v). Otherwise, is the current 

number of unsatisfied constraints that can be satisfied if the free 
variable j is fixed at 1. Initially, n(0,j) = ^ a., for all 

j. ■'=1 

The current number of unsatisfied constraints is kept in a vector <}>, i.e., 
<f)(o) = m and (j)(£) - 4>{ii'-l) - ri(Jl - 1, j*) . For any £,j, the condition 
<j)(Jl) > n(.'J-,j) ^0 is always held. 
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Consider the search at level JZ- prior to a forward step. The following tests 
may reduce the number of branch alternatives. 

n 


a. The subproblem at level a is feasible only if 
this is not met, a backward step may be taken. 


j=l 


b. Let Z be the current objective function (2^ °^er all j with y^ 
fixed at 1) and Z* the objective function of the best integer solution. 
Let F be the set of free variables. Then the double ceiling test, i.e., 

^ pennits the cancelation of t if 

n(£,t) < 4)(s,). 


c. If at a current subproblem, a variable column dominates another variable 
and the dominated column has a higher cost than the dominating one, then 
the dominated variable never need be considered as a branch candidate, 
that is, if for any level Jl{jl^l), n(Jl - 1, j*) ^ n(A - 1, j-j) > 0 and 
= n(iZ-,j-j) = 0 with then set n(Jl - 1 , j-j) = 0 {i.e,, 

cancel at level £ - 1). Here y.* is the branch from node v" to 
V. The above is commonly referred to as the local column dominance test. 


^ a 

d. Let (I) be the level Jl Integer (set covering) subproblem. Define 
(Ilr)^ to be the corresponding continuous LP problem. Then as stated 
in section 2.3.1 ; 


{i) feasible if and only if (I)^ is feasible. 

(ii) If (I, n)^ is feasible, then the value of the minimal functional 

Ur .s 

serves as a lower bound for the best (I) solution. Furthermore, 
if (I|_p)^ is integer feasible it solves (I)^. 


2.3.3 THE ALGORITHM 

The algorithm is basically a (0,1) search embodying the elaborations discuss- 
ed previously. At node v, level I, we attempt to solve the linear program 
(Ij_p)®'. If it is infeasible or its objective function exceeds the current 
ceiling (z*, initially set to c'e) the search reverts to v*. If the minimal 
solution is (0,1) feasible, it is recorded and a backward step is taken. 
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Otherwise, we record the column (denoted by j*) associated with the smallest 
value of the optimal variables, and then extract the sub- subproblem that 
contains the columns of A associated with the positive fractional variables 
and the rows corresponding to the constraints they satisfy, and obtain a 
solution to this problem. This solution should be extreme point which is 
obtained through the procedure introduced in 2.3.1. If the overall (I)- 
feasible solution exceeds the ceiling z*, a backward step is taken. 


If node V was just reached on a forward step, we construct the preferred 
set P{i*). Preferred row i* is defined to be the first i such that 

SjeF ^i*j -^jeF ^ij which j*th entry is 1. P(i*) is 
the columns containing 1 on the preferred row i*. We then cancel locally 
dominated columns. The purpose of constructing the preferred set is if a 
forward step fails, we cancel j* and select another branch from the pre- 
ferred set. 


If V was reached from v', we branch j* to reach node V'5', However, if 
V was reached on a backward step (that implies when we tried to reach 
and found that we cannot generate an improved solution), we select a branch 
by the minimum-cost-per-constraint-satisfied rule from P(i*) to define v"^. 
During the branch selection procedure, the double ceiling test is applied so 
we can hofefully exclude some branch candidate from further consideration. 

If at any time every variable in the P(i*) is canceled, we take a backward 
step. 

When a subprogram contains no constraints, we have reached an (I)-feasible 
solution. In all cases, a (0,1) solution is not recorded until it is reduced 
to an extreme point of (iLp)* The search terminates when the level counter It 
becomes less than 0. 

Linear programming is the crucial feature of this algorithm. Initially, the 
sub-subprogram is always extracted if (I^p) Is not integral feasible. A 
solution of the sub-subproblem is reached by successively setting variables 
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to 1 by the "minimum-cost-per-constraint-satlsfled" rule. Once this solution 
is reached, we branch forward by setting the minimal fractional variable of 
the current LP to 1. The hope is that setting the minimum fractional vari- 
ble to 1 will tend to give the largest alteration in the optimal LP solution 
which is the lower bound of (I). When a node is reached for the first time 
on a backward step, LP is performed with the minimum fractional value at 0. 
This means that last branch link j* has been canceled. 

The other tests, such as local dominance, double ceiling, and so on, were 
incorporated more because of their simplicity and ease of applicability than 
their actual usefulness in curtailing the length of the search. In dealing 
with the usual large number of variables, search features designed to look at 
individual variables appear relatively less than attractive as compared with 
the linear programming part that tends to solve entire subproblems. This is 
reflected by the fact that for each computer run, 90 percent of the computing 
time is consumed in LP, The saving of search time because of successful 
LP solution appears more important than the tests mentioned above; especial- 

/N 

ly after we experienced most of the time that the (I)-feasible solution was 
directly generated by LP. 

A general functional flowchart is given in section 5.1. 

2.3.4 OUTLINE OF THE ALGORITHM 
I. Initialization 

1.1 Set n{l»j) = i"or j = 1, ...,n; F = {l,2,...,n}; s = e 

1.2 = m,?(l) = -1, o{l) = 0 

1.3 J!, = 0, 7 = 0, z = 0, z* = c*e, a = 1 

II. Point Algorithm 

2.1 Updating counters + = 

2.2 Test for solution and lower bounds. If 4»{A) == 0, go to 4.1 

(solution). If ^ z* -T, go to 3 (backward step). 

2.3 Update 0, F. Set + 1) = (!)(£.), F - {j i ) |>0}. (Variables 
with n(A,j) entries less than 0 are r.\embers of the preferred set 
See 2.12). 
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2.4 Feasibility test. If |ti(S-,o) |«t>(5') » go to 3. 

i=l 

2.5 Linear programming. Solve the linear program associated with the 
current subproblem. 

(i) If infeasible or ceiling exceeded, go to 3. 

(ii) If (0,1) feasible, record overall (I) solution, go to 4,1. 

(iii) Otherwise, set ?(il) equal to the minimal LP objective function 
added to the current value of z; and set j* corresponding to the 
smallest fractional value of the LP variables. 

Obtain the roundup solution if currently feasible, and go Lo 2.6; 
otherwise, go to 2.7 when a ~ 0 and go to 2.10 if a = 1. 

2.6 Sub- subproblem. Extract and obtain a solution to the sub-subproblem. 
Record the overall (1) solution. Go to 4.2. 

2.7 Feasibility test. If any j I = 0. set j* = j (solution 

at next level), and go to 2.13. 

2.8 Double-ceiling test. If any j such that n(Jt.,j) > 0, 

2 + Cj + i " '^9 set n(5-,j) " 0 (i.e, , cancel y. at the 

level if n()l,j) = 0 for all j, go to 3, 

2.9 Select branch j(a = 0), If ’n(fi.,j) > 0 for all j (preferred set 
null), go to 3. Set j* = first j to satisfy 

l] < [':j/[n(^,j) I j for all j such that n(fi',j) < 0 
(i.e,, for all j in the preferred set). 

2.10 Update n. Set cj).(£ + 1 ) = + 1 ) - ln(i!..j*) I , 

n(S- + Ij) = !n(^l,j)| for j = l,..,n. If (j)(a + 1) = 0, go to 2.13. 
For each i such that < 0 (i.e., for all rows of the subproblem) 

and a..* = 1, set n(J^ + U) to n(2- + l>j) - 1 when a.. = 1 

I * J 

and jeF. 

2.11 Local dominance test. If, for any jeF, n(2- + 1>J) = 0 and 
Cj* < c j , set ri{)l,j) = 0, If a = 0, go to 2.13. 

2.12 Construct preferred set. Set (the preferred rovj) to be the first 

i such that Sjgp - ^jcF ®ij ^ ^i ^ 

n(Jl,j) to -n(J?-,j) for the j such that a.^^. = 1. 

2.13 Update the slack column and other parameters. Set s^ to s^. + 
for i = l,...,m, z = z + c.^, a = 1, ?(.?, + 1) = C(Z). Go to 2,1. 

tJ 
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III. Backward Step 

3. Set l-l ~ Z. If it < 0, go to 5 (termination), j* = cr(il + 2). 

z = z - c.*, s. - s. - a..* for i = + 2) = 1). 

a = 0. Go to 2.1. 

IV. Feasible Solution 

4.1 LP (0,1) feasible, or search reaches solution. Reduce overall 

solution to extreme point. Record current solution y and 
c'y = z. Set z* = z. Go to 3. 

4.2 Search produced via sub-subproblem. Reduce overall solution to an 
extreme point. Record the overall y and c'y = z. Set z* = z. 

If C(iJ') > z* - r, go to 3. If a = 1 go to 2.10, otherwise 

(a - 0) go to 2.7. 

V. Termination 

5. Optimal solution ascertained or no feasible solution exists. 

2.3.5 AN EXAMPLE 

The following example will demonstrate how the algorithm is applied to the 
traffic model problem. The input obtained from MPLS is an occurrence table 
which is as shown in table I. 

The constraint matrix is displayed as in table II. Assume we intend to find 
a traffic model with the minimum number of missions and the cost coefficients 
are in the column matrix of all I's, Since the traffic model is a set 
partition problem, we transform the cost function as indicated In section 
2.3.1 with M = 100. The new cost column is as shown in table II. The 
search follows the outline of the Set Covering Algorithm. As an example, the 
enumeration of the problem in table IT is demonstrated in table III. In each 
level on the forward step, the LP subproblem is solved and the lower bound 
is updated. Ceiling z* is recorded after the rounded-up solution of LP 
has been reduced to an extreme point. The search always branches forward 
until it reaches node 6. At node 6, the lower bound is equal to the ceiling, 
so a backward step is taken. In each backward step, it is found that 
z* - < T, so z* is the solution. 
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TABLE I.- OCCURRENCE TABLE 


Payloads 
1 
2 
. 3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


Feasible combinations 
1 , 16 , 17 , 18 . 19 , 20,21 , 64 , 65 , 66 , 67 , 68 , 73,74 

2 . 16 . 23 . 24 . 25 . 26 . 27 . 28 . 29 . 64 . 65 . 66 . 69 . 70 . 71 . 73.74 

3 . 17 . 23 . 30 . 31 . 32 . 33 . 34 . 35 . 64 . 67 . 68 . 69 . 70 . 73.74 
4,36 

5 , 37 , 38 , 39 , 40,41 

6 , 37 , 42 , 43,44 

7 , 38 , 45 , 46,47 

8 , 48 , 49 , 50,51 

9 , 48 , 52 , 53,54 

10 , 22 , 24 , 30 , 55,71 

11 . 18 . 25 . 31 . 39 . 42 . 45 . 49 . 52 . 56 . 57 . 58.72 
12 , 19 , 26 , 32 , 59 , 60 , 65 , 67 , 69,73 

13 . 20 . 27 . 33 . 36 . 55 . 56 . 61 . 62 . 66 . 68 . 70 . 71 .74 
14,21 , 28 , 34 , 40 , 43 , 46 , 50 , 53 , 57 , 59,61 , 63,72 

15 . 22 . 29 . 35 . 41 . 44 . 47 . 51 . 54 . 58 . 60 . 62 . 63.72 
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TABLE II.- AH EXAMPLE CONSTRAINT MATRIX 


jA blank entry In the constraint matrix means zero] 


Col. 1 2 3 4 6 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 

c' (101 TOT 10! 101 101 101 101 101 101 101 101 101 101 101 101 201 201 201 201 201 201 201 201 201 201) 


Row 







TABLE II.- AN EXAf^PLE CONSTRAINT MATRIX (Continued) 
[A blank entry in the constraint matrix means zero] 


Col. 

Z6 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

e' 

(201 

201 

201 

201 

201 

201 

201 

20T 

201 

201 

201 

201 

201 

201 

201 

20i 

201 

201 

201 

201 

201 

201 

201 

201 

201} 


Row 







TABLE n.- AN EXAMPLE CONSTRAINT MATRIX (Cont'inued) 
[A blank entry In tbe constraint inatrix means zero] 


Col. 51 52 S3 54 55 5fi 57 58 59 60 61 62 63 64 65 66 67 68 69 7 0 71 72 7 3 74 

c' (201 201 201 201 201 201 201 201 201 201 201 201 20l 301 301 30V 301 301 301 301 301 301 401 401) 










































































2.3.6 POSSIBLE IMPROVEMENTS 

The following modificati'^ns may result in the improvement of program 
efficiency. 

1, Based on the fact that the rounded-up solution is always feasible, in 
outline 2.5, after obtaining the rounded-up solution, control should 
always go to 2.6. 

2. In the outline 2.11, the local dominance test may lead to an infeasible 

I 

solution unless we stress column j is dominated by column j*. For 
example, consider the problem 

min 5x, + X 2 + 2x^ 


X, + Xj X3 =1 

Cg = 1 < Cl = 5 and assuming j* = canceling x^ will lead to an 
infeasible solution which will not satisfy the firsL v-ow of constraints. 

3. In solving a set partition problem, any jeF such th.-'t A. • aT* 0 

w vJ 

should be canceled on a forward step. By doing that, not only can we 
exclude more branch candidates, but also we can save the computing time 
in solving the LP subproblem. For example, in tables II and III, in 
the first level according to the algorithm we only canceled y.^^ and 
y.jg. However, noticing the equality constraints, yg^, y^g, y^g, y^i* 

^29* -^35* ^4r -^44’ ^47, ^54’ ^58’ ^60* ^62’ ^63 ^72 

canceled also, since any of those variables not equal to zero will result 

in an infeasible solution by branching on y 22 - 
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3. PROGRAM USAGE 


3,1 INPUT DESCRIPTION 

The two types of input data for the SCA program are the source input from a 
file on logical unit 2 and the tutorial input data specified by the user. 

All source input data are nonnegative integars. Logical unit 2 is built as a 
temporary file in the main program to store the feasible payload combinations 
from the MPLS. Because of the core storage limitations, the SCA is designed 
to handle a traffic model problem with a maximum of 500 combinations and 100 
payloads. If the number of feasible combinations exceeds 500, the MPLS will 
reduce it to within that limit. Subroutine TABLE reads in the information on 
each feasible combination from logical unit 2. Within the information combi- 
nation ID, different payloads in that combination, and combinations cost co- 
efficients are stored. All the inputs the SCA needs is a vector c, the vec- 
tor of cost coefficients, and a constraint matrix A which is displayed as 
an occurrence table. The cost vector is directly stored in the column array 
KCJ. The matrix A is stored as two vectors, i.e,, KARR and NOPERC. For 
example, suppose A is 

column 1234 
row 
1 
2 
3 

Then A is posed as a vector KARR q ~ (1 >3;2;1 ;2,3), which contains the 

indices of rows of all nonzero entries columnwise from left to right. The 
position vector, NOPERC (I)j=i 5 = (1,3, 4, 5, 7) containing n + 1 elements, 
keeps track of the column interval in KARR. This implies that column j of 
A starts in KARR (NOPERC(j)) and ends in KARR (N0PERC(d+l ) - 1). For 
example, we want to locate the Vs in column 4. 

Set j = 4, then KARR(NOPERC(j) ) = KARR(5) =2, 
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KARR(N0PERC(j+1) - 1) = KARR{7-1) = KARR(6) = 3; so column 4 hcs two Ts, 
on row 2 and 3, respectively. 

The tutorial input data required for the 5CA program can be either user- 
specified data from a demand terminal or from card decks. These data are 
read via logical unit 5 using a free field format. Sample input/output is 
given in section 5.4 for the reader's reference, 

3.2 PROGRAM RUN PREPARATION 

The SCA has been implemented on the UNIVAC 1110 EXEC 8 system as a subprogram 
of SAMPLE. FORTRAN V standard logical input and output devices are used for 
tutorial input (logical 5) and printed output (logical 6). For the source 
input of the SCA, logical 2 is used. 

In order to eliminate the reiteration of usage instructions which have already 
been elaborated in the SAMPLE User's Guide (ref. 2), a discussion of the 
natures of those different interactive options will be introduced as follows 
instead. 

3.2,1 INTERACTIVE OPTIONS 

1 : USE PREVIOUSLY DEFINED FEASIBLE COMBINATIONS 

2: USE INTERACTIVE FEATURE IN TRAFFIC MODEL 

3: NONE OF THE ABOVE 

Option 1 will enable the user to use a data file input which contains previ- 
ously defined feasible combinations for SCA execution. In this way, the 
execution time of the MPLS can be saved. 

Option 2 mainly supports the communication between the SCA and the OSARS. In 
case the traffic model cannot be scheduled by the OSARS, or is not desired 
for some other reason, this option gives the user the means of changing the 
traffic model. 

Option 3 implies that the user does not desire to select option 1 or 2. 
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3.2.2 COST CRITERIA OPTIONS 

This option gives user the choice of one of the performance criteria against 
which traffic model will be generated. The meanings of different criteria 
have been discussed in section 2.2,2. 

3.2.3 CRITERIA FOR FLIGHT/ COMB I NAT I ON SELECTION OPTION 
CHOOSE CRITERIA FOR FLIGHT/COMBINATION SELECTION: 

j 

1 : MAXIMUM NUMBER OF PAYLOADS 

2: MAXIMUM PRIORITY 

3: MINIMUM COST 

4; MINIMUM COST PER PAYLOAD 
5: NONE OF THE ABOVE 

These criteria help the user to make the decision about which mission he likes 
to enter in the traffic model. If the user is interested in adding certain 
feasible combinations which contain the largest number of payloads* he can 
select option 1. Then the program will print out five feasible combinations 
with maximum number of payloads. Presently, the priority is determined by 
the number of payloads, so option 1 and option 2 give the same output. By 
selecting option 3, the user will get a list of five feasible combinations 
with minimum cost coefficient. Option 4 will give a list of five feasible 
combinations cost coefficient per payload. Control goes to the Manual Flight/ 
Combination Option when 5 is chosen. 

3.2.4 MANUAL FLIGHT/COMBINATION OPTIONS 

The available options are: 

N: ENTER COMBINATION "N" 

-1 : NEW SELECTION CRITERIA ARE DESIRED 

-2: CONTINUE ON TO TERMINATION 

-3: VIEW ALL COMBINATIONS SPECIFIED SO FAR 

-4: VIEW INFORMATION ON COMBINATIONS SPECIFIED SO FAR 

-5: REMOVE LAST SPECIFIED COMBINATION 
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This part follows right after "Criteria for Flight/Combination Selection 
Option" in which the user has viewed the relevant information on the combi- 
nation he possibly adds in the traffic model. If the user would like to 
add a specific mission in the traffic model, he responds by entering the 
feasible combination number "N." The program will take that feasible combi- 
nation into the SCA's partial solution. If the user still wants to see more 
relevant criterion selection information, he just enters -1. The control 
goes back to Criteria for Flight/Combination Selection Option. Option -2 
makes the control continue to find a traffic model with the specified 
partial solution. A traffic model is built around the partial solution; 
local optimality replaces global optimality. Option -3 gives the user a 
chance to lo.ok through all the missions already specified in the partial solu 
tion so he will not enter any of those missions again. By selecting option 
-4, the user will see a list of missions specified in the partial solution 
and their relevant information on Shuttle sequence, inclination, payload 
margin, and so on. After the user looks through the detailed information on 
those specified combinations and he is not satisfied with the mission he 
just added on the partial solution, he can enter option -5 to remove it. 

3.2.5 TRAFFIC MODEL INFORMATION OPTION 

The tutorial of this option is 

DO YOU WISH TO SEE INFORMATION ON THESE MISSIONS? 

0: NONE 

-1 : PRINT ALL 

-2: PRINT ALL AND SAVE ON SCRATCH FILE 

-3: SAVE ON SCRATCH FILE ONLY 

N: ENTER MISSION "N" • 

This option follows the solution of the traffic model. If the user does not 
need to see any detail information on the traffic model, he just enters 0. 

He can enter -1 to see them all or enter N, the combination number, to get 
the information on a particular one. Options -2 and -3 provide the user an 
opportunity to store the traffic moJel on a scratch file (logic unit 1) for 
further analysis. 
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3.2.6 MISSION OMIT OPTION 

The display statement is 

WHICH MISSIONS DO YOU WANT OMITTED? 

After the user viewed some information on the traffic model, he may need to 

delete certain combinations by use of this option. 

3.2.7 TERMINATE OPTION 

The display statement is 

SELECT AN OPTION: {3 TO TERMINATE) 

Input 3 to terminate execution. 

3.3 OUTPUT DESCRIPTION 

3.3.1 NORMAL OUTPUT 

Normal output for the SCA program can be classified into five basic types: 

1. Source input data - The initial output of the SCA is the source input 
data which is displayed in the occurrence table. The title of that table 
is printed out as "n OCCURRENCE TABLE," where n is the year with which 
the particular case is executed. This is immediately followed by 
"PAYLOAD" and "COMBINATIONS." Under the column of "PAYLOAD" are printed 
out the payload identifications. Under the column of "COMBINATIONS" are 
combination numbers which carry that payload. 

2. Tutorial instructions data - These data are printed out in the alpha- 
numeric format and provide the user a guide of various interactive se- 
lections during the execution of the SCA. 

3- Criteria for flight/combination selection data - The output is written 
out in two columns; the first column contains the payload identification, 
and the second column displays the corresponding criterion. 

4. Combination information data - The output of these data is in alpha- 
numerical format and displays the relevant information about the mission 
in an understandable form. This output is requested by the user In the 
Manual Flight/Combination Option. 
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5. Traffic model data - These data give the total number of combinations in 
the traffic model, the mission identifications, and the total cost of 
this traffic model. 


All of the five output types can be identified in a sample input/output in 
section 5.4> 

3.3.2 ABNORMAL OUTPUT 

Diagnostic messages from subroutines of the SCA are listed below. 


Diagnostic Message Subroutine 

THIS MISSION IS UNACCEPTABLE SET 

BECAUSE PAYLOAD XX IS DUPLI- 
CATED 


ERROR, ALL ROWS ARE NOT SET 

COVERED, RETURN TO PREVIOUS 
SOLUTION, PLS NOT COVERED 
ARE XX, YY, ZZ 


REDIMENSION KARR (XX) TO SET 

KARR (YY) 


Description/Action 

In the flight/combi- 
nation option, the user 
entered more than one 
combination which covers 
the same payload. XX is 
the payload ID, 

The user wanted to omit 
some combinations that 
will cause some payloads 
not to be covered. The 
program will return to a 
previous partial solu- 
tion. If the user wants 
to omit some combinations 
this time, he should 
refer to the occurrence 
table and be sure all 
payloads can be covered, 
XX, YY, ZZ are the pay- 
loads which failed to be 
covered. 

The total number of unity 
entries in the con- 
traint matrix YY which 
exceeds the dimensioned 
space XX. 
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Diagnostic Message 


Subroutine 


Description/Action 


AVAILABLE STORAGE EXCEEDED SET 

AFTER ITERATIVE STEP NO. I 


TABLE ERROR***INPUT TO TABLE 

SCA IS CLOBBERED** 


Variable KX8 (I,J) 
should be redimensioned 
by increasing the value 
of 1. The value of 
KTEST should also be 
increased by the same 
amount. 

This message impne^ 
that more than one pay- 
load in a feasible 
combination have the 
same ID. It could be 
caused by the numbering 
or naming method in the 
MPLS. 
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4. EXECUTION CHARACTERISTICS 


4.1 RESTRICTIONS 

The SCA program has these limitations: 

a. The largest traffic model problem the SCA can accept is 100 payloads with 
500 combinations. 

b. The level the SCA can reach is limited to a maximum of 6. 

c. The program is valid only if there exists a feasible solution to the 
traffic model, 

d. The maximum number of iterations allowed 1n LP is 243. 

4.2 RUNNING TIME 

The run time for the SCA program may vary depending on the problem executed. 

A rough estimate of the time needed for a run can be obtained from the plot 
of the number of missions as a function of running time in figure 2. The 
data from which this plot is constructed are from 12 cases with unity cost 
coefficients. 

4.3 ACCURACY/VALIDITY 

The SCA program is written in single precision, and has been checked out 
using small problems, in a range from 11 feasible combinations and 3 pay- 
loads to 229 combinations and 97 payloads, with a maximum of 3 payloads 
per combination. It is felt that the program is operating correctly and 
is providing reliable solutions to the problems. 

Problems used in checkout included three data sets and 36 problems. Each 
of these problems has been tested on four different sets of cost coeffi- 
cients with a maximum of three payloads per combination. All solutions 
have been checked and were found to satisfy the constraints and to be 
optimal. Some other problems with more than 400 combinations which have a 
maximum of four payloads per combination were also run and resulted in an 
indication that no solutions were possible. The reason for this has not 
been thoroughly analyzed. 
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T^UN TIME IN SECONDS 


ao l*Maximum number of payloads per combination - 3 


70 I 


A 


60 L 


50 L 


40 


30 {■ 


20 L 




A 


A 


10 ^ 



0 5C 100 150 200 250 300 

NUMBER OF FEASIBLE COriBINATIOHS 
Figure 1.- Sample run times. 
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5. REFERENCE INFORMATION 


5.1 FUNCTIONAL FLOWCHART 

Figure 2 illustrates the flow of the controlling subprogram logic. Refer 
to section 2.1 for the definitions of those symbols used in this functional 
flowchart . 




Figure 2.- SCA functional flowchart. 
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Set equal to 2 + (minimal LP 
object function). Set next branch j* 
corresponding to the smallest frac- 
tional value of LP variables. Obtain 
a round up solution. Extract the norv| 
integer variables and associated 
constraints to form a sub-subproblem 


In subroutine PURIFY, obtain a solution of| 
the sub-subproblem. Record the overall 
integer solution. Reduce overall solution 
to an extreme point. Record the solution 
y and c'y = z. Set z* = 2. 


Yes 





Figure 2.- Continued. 
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Figure 2.- Continued. 
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Figure 2.- Continued, 
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Figure 2 Concluded. 
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5.2 SYMBOL DEFINITIONS 

Table IV defines a list of parameters specified in DATA statements in the 
SCA subprogram. Table V gives the description of all variables used in 
labeled COMMON. 
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TABLE IV - PARAMETERS IN DATA STATEMENTS 


Parameter 

name Dimension Type Value 

KSTRl 1 I 1 (0) 


KSTR2 1 I 1 (0) 


KSTR3 1 I 1(0) 


KSTR4 1 I 1(0) 


Description 

When KSTRl = U the 
algorithm always extracts 
and gets a (0*1) solution 
from the sub-subproblem. 
Otherwise, the sub-subproblem 
is explored only when the 
LP roundup solution is 
currently feasible. 

When KSTR2 = 1 , the 
algorithm will attempt to 
reduce (0,1) feasible 
solutions to (0,1) feasible 
extreme points of the 
associated continuous 
problem. Otherwise, this 
attempt will not be made. 

When KSTR3 = 1 , the search 
will select branches by 
the minimum-cost-per- 
constraint-satisfied rule; 
otherwise it uses the 
maximum-number-of-constraints- 
satisfied rule. 

When KSTR4 = 1 , the search 
will ignore other criteria 
and branch on the minimum 
fractional LP value when 
available- Otherwise, the 
branch selection rule is 
determined by KSTR3. 
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TABLE IV. - PARAMETERS IN DATA STATEMENTS - Concluded 


Parameter 

name 

Dimension 

Type 

Value 

Description 

KSTR5 

1 

I 

1 

When KSTR5 = 1 , 


algorithm supposes that 
the user will supply a 
slack cost, either by 
setting it directly to MM 
or by supplying MAXCST, in 
which case MM is set to 5* 
MAXCST. If KSTR5 = 0, the 
algorithm automatically 
takes the slack cost as 0. 

KSTR6 1 I 1 When KSTR6 = 1 , LP is 

performed after forv/ard 
and backward steps. Other- 
wise LP is performed only 
after forward steps. 
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TABLE V. “ VARIABLES IN LABELED COMMON 


« COMMON Block Name: C9 

Description: C9 retains the information about the interactive selection of 
the cost coefficient of the objective function and the output of the occur- 
rence table. 


Location * 

Name 

Dimension 

Jm. 

Description 

2 

MM 

1 

1 

Total number of feasible combi- 
nations generated by MPLS 

51 

COSTOP 

1 

I 

Indicator of the choice of cost 
coefficients on the objective 
function 

54 

NOTAB 

1 

I 

Indicator of listing or suppressing 
the occurrence table output 


• COMMON Block Name: Cl 3 

Description: C13 retains information about the structure of the constraint 

matrix which defines the traffic model problem arsd a flag to trigger the OSARS. 


Location 

Name 

Dimension 

Type 

Description 

1 - 501 

NOPERC 

501 

I . 

The position vector indicating the 
number of ones in each column of 
the constraint matrix 

502 - 2501 

KARR 

2000 

I 

The position vector indicating the 
rows which correspond to one entries 
in the constraint matrix 

2502 

KM 

1 

I 

The number of rows of the constraint 
matrix 


*The unspecified cells are used for the OSARS. 



TABLE V. - VARIABLES IN LABELED COMHON - Continued 


• COMMON Block Name: 

Cl 3 - Continued 


Location 

Name 

Dimension Type 

Description 

2503 

KN 

1 I 

The number of columns of the 
constraint matrix 

2504 

NOSARS 

1 I 

A flag indicating the user’s 
choice of the use of the OSARS 


# COMMON Block Name: Cl 6 

Description: Cl 6 contains information about the location of testing range of 

each combination. 

Location Name Dimension Type Description 

1 - 500 LOCEOW 500 A A vector to store alphanumerical 

identification of testing site 
where a particular mission is to be 
launched 


9 COMMON Block Name: CIS 

Description: CIS contains information about the coefficients and the priorities 

of those variables in the objective function. 


Location 

Name 

Dimension 

Type 

Description 

1 - 500 

KCJ 

500 

I 

The cost coefficients of the ob- 
jective function is to be mini- 
mized in the SCA. They may be 
defined as feasible combination's 
OMS weight, one-load factor, or 
one- payload length 

501 - 1000 

KPRIOR 

500 

i 

The weight factor to determine 
each feasible combination's pri- 
ority. 



TABLE V.- VARIABLE IN LABELED COMMON - Continued 
e COMMON Block Name: C27 

Description: C27 retains various information about the branch-and-bound 
algorithm used in the SCA. 


Location 

Name 

Dimension 

Type 

Description 

1 - 50 

ZLB 

50 

F 

The lower bound of the 
branch-and-bound method 

51 - 3550 

KXB 

7 X 500 

I 

The number of unsatisfied 
constraints of each variable 
at different levels 

3551 - 3650 

KYS 

100 

I 

The value of the slack variabl 

3651 - 4150 

KCOM 

500 

I 

A vector to store the column 
number of those variables in 
the preferred set 

4151 - 4651 

KPREF 

501 

I 

A vector to store the row 
number of the infeasible 
constraints 

4652 - 5151 

LBSC 

500 

I 

An indicator of a decision 
variable's value in the LP 
solution. LBSC(K) negative 
means variable K is one. 

5152 - 5651 

LP3TA 

500 

I 

A vector to store the pivot 
row of each LP iteration 

5652 - 6252 

LI DRW 

601 

I 

Scaled cost coefficients of 


the objective function 
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TABLE V.- VARIABLE IN LABELED COMMON - Continued 


• COMMON Block Name: C28 

Description: C28 transmits information between the 

SIMPLX routine. 


Location 

Name 

Dimension 

Im. 

1 

KPV 

1 

I 

2 

KHP 

1 

I 

3 

KTEST 

1 

I 

4 

MM 

1 

I 

5 

NLP 

1 

I 

6 

KSOL 

1 

I 

7 

KUNF 

1 

I 

8 

KFEAS 

1 

I 

9 

KZS 

1 

I 

10 

KZSTR 

1 

I 

11 

J2 

1 

I 


SET routine and the 

Desc ription 

Indicator of the number of 
elements in the preferred set 

Level indicator of the SCA 

The highest level allowed in the 
SCA 

Cost coefficients of the slack 
variable of the objective function 

An indicator of the number of 
times subroutine SIMPLX has been 
called 

Indicator of the number of 
improved solutions reached in the 
SCA 

A variable to indicate the feasi- 
bility of the LP solution 

An indicator of the binary feasi- 
bility of LP solution 

Current value of the objective 
function 

Current upper bound of the 
objective function 

Branching variable in the SCA 
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TABLE V.- VARIABLE IN LABELED COMMON - Continued 
« COMMON Block Name: C28 ^ Continued 


Location 

Name 

Dimension 

Jm. 


12 

K1 

1 

I 

The number of variables in the 
preferred set 

13 

NEP 

1 

I 

The number of times an extreme 
point reduction option is used 

14 

NEPSC 

1 

I 

The number of times an extreme 
point reduction is successful 

15 

NOGO 

1 

I 

A flag indicating infeasibility 
in subroutine PURIFY 

16 

LD7 

1 

I 

A flag to print out final inf or 
mation for subroutine PURIFY 

17 

KSUM 

1 

I 

Current value of the objective 
function in a sub-subproblem 
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TABLE V.- VARIABLE IN LABELED COmON » Concluded 


ir COMMON Block Name; C29 

Description: C29 retains various inform?<;t1on about the partial solution of 

the SCA. 


Location 

Name 

Dimension 

Im. 

Description 

1 - 500 

KSEQU 

500 

I 

A vector to store the partial 
solution 

501 

KSC 

1 

I 

Level indicator 

502 

KANDO 

1 

I 

A flag indicating that the OSARS 
can schedule the traffic model when 
it equals one 

503 

KNSKDl 

1 

I 

A flag to choose the solution 
output format for the SCA or 
the OSARS 

504 - 553 

KPSKED 

50 

I 

A vector to store the missions 
scheduled by the OSARS 

554 

KPRMAX 

1 

I 

Best attainable mission priority 

555 

KNOSKD 

1 

I 

The number of missions scheduled 


by the OSARS 
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5.3 SUBROUTINE DOCUMENTATION 


Individual subroutine documentation appears in alphabetical order on the 
following pages. 
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SUBROUTINE BISRCH 


IDENTIFICATION 

Name/Title 
Author/Date 
Machine Identification 
Source Language 


-- BISRCH (Binary Column Search) 

- Han Chang, July 1975 

- UNIVAC 1110 

- FORTRAN V 


PURPOSE 

Subroutine BISRCH searches and positions particular row entries for the 
given column of a constraint matrix. It was written to update the n 
matrix, generate a preferred set, test for an extreme point, etc. 

USAG E ' 

« CALLING SEQUENCE 
CALL BISRCH C$,K,J) 


Arguments: 


Parameter 

Name 

In/Out 

Dimension 

Iffie 

Description 

$ 

Out 

1 

I 

Nonstandard return signal 
when a particular entry 
has been found. 

K 

In 

1 

I 

The column index on which 
a particular row entry is 
to be searched 

J 

In 

1 

I 

Row index on which an entry 
is to be positioned at 
a given column K 


BISRCH-1 


5/75 
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fi Data In/Out 

Labeled COMMON (refer to the labeled COMMON description section): 


Block Name 

Input 

Output 

CIS 

1 - 50T 
502 - 2501 



METHOD 
e Model 

Subroutine BISRCH searches over a column interval of the constraint matrix 
to locate a -particular entry. The output consists of a nonstandard return 
whenever this entry has been found. 


5/75 


BISRCH-2 
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SUBROUTINE FNDFLT 


IDENTIFICATION 

Name/Title 
Author/Date 
Machine Identification 
Source Language 

PURPOSE 

Subroutine FNDFLT provides the necessary information about a combination 
by the interactive request from the user. It also saves that information 
from the traffic model in a scratch file of the user’s choice. 


- FNDFLT (Find Flight) 

- Han Chang, July 1975 

- UNI VAC 1110 

- FORTRAN V 


USAGE 

9 CALLING SEQUENCE 
CALL FNDFLT (KZS,LIDRW,IPONT) 

Arguments : 

Parameter 


Name 

In/Out 

Dimension 

Type, 

Description 

Total number of missions to 

KZS 

In 

1 

I 

be output 

LI DRW 

In 

1 

I 

A vector to store the mission 
numbers 

I PONT 

In 

1 

I 

A flag to trigger the output 


of the statistics of the 
current flight schedule 


8/75 


FNDFLT-1 
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s DATA In/Out 

Labeled COMMON (refer to the labeled COMMON description section): 

Block Name Input Output 

Cl 3 2503 

METHOD 

Subroutine FNDFLT searches over a data file (logic unit 2) to locate parti- 
cular missions specified by the user. Detailed information on those combi- 
nations are output in an understandable format by calling subroutine DISPLY. 
Information about the combinations in the traffic model is saved on the 
scratch file (logic unit 1) at the user's request. 


8/75 
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SUBROUTINE PURIFY 


IDENTIFICATION 

Name/ Title 
Author/Date 
Machine Identification 
Source Language 


- PURIFY (Purification of LP Solution) 

- Han Chang, July 1975 

- UNIVAC 1110 

- FORTRAN V 


PURPOSE 

Subroutine PURIFY purifies noninteger linear programming solution into an 
improved rounded-up integer solution. 


USAGE 

e CALLING SEQUENCE 
CALL PURIFY 

e Data In/Out 

Labeled COMMON (refer to the labeled COMMON description section): 


Block Name Input Output 


Cl 3 

1 

- 501 


502 - 

2501 



2502 

CIS 

1 

- 500 

C26 


4 

C27 

3551 - 

3650 


4151 - 

4652 


4653 - 

5154 


5155 - 

5654 


5655 - 

6255 


5/75 
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Block Name 


Output 


Input 

C28 , 1 

2 

6 

9 

10 

13 

15 

17 

C29 1 - 500 

501 


15 

17 

1 - 500 
501 


METHOD 
e Model 

Wtienever the linear programming solution is not integrally feasible, a sub- 
problem is created by eliminating the columns associated with variables 
having an integer value and those constraints satisfied by these variables. 

In the subproblem, PURIFY successively sets to one the variable corresponding 
to the minimum'*cost-per-constraints-satisfied ratio until all constraints 
have been satisfied. 

REFERENCE 

C. E. Lemke, H. M. Salkin, and K. Spielberg, "Set Covering by Single 
Branch Enumeration with Linear Programming Subproblems," Operations 
Research 19, pp. 998-1022 (1971). 
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PURIFY-2 

5-22 


SUBROUTINE RIDMOD 


IDENTIFICATION 

Name/Title - RIDMOO (Read Input and Mod) 

Author/Date - Han Chang, July 1975 

Machine Identification - UNIVAC 1110 

Source Language - FORTRAN V 

PURPOSE 

Subroutine RIDMOD provides the selection of flight/combination information 
for user's reference when he goes through the SCA interactively. Based on 
this information, the user can make his judgment as to which combination 
he would add or delete from the traffic model. 

USAGE 

« CALLING SEQUENCE 
CALL RIDMOD 

« DATA In/ Out 

Labeled COMMON (refer to the labeled COMMON description section): 

Block Name Input Output 

C31 1-4 


METHOD 

Subroutine RIDMOD collects the user's numerical choices on flight/combination 
option and sorts them out by using a MOD fu’ "tion. Then subroutine SET uses 
this information as input to generate the p, Jp'r flight/combination information. 


5/75 


RIDMOD-1 
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SUBROUTINE SET 


IDENTIFICATION 

Name/Title - SET (Set Covering Algorithm) 

Author/Date - Han Chang, July 1975 

Machine Identification - UNIVAC 1110 

Source Language - FORTRAN V 

PURPOSE 

Subroutine SET is the main driver of the SCA. SET determines the feasibility 
of the intermediate solution from subroutine SIMPLX on each level and decides 
whether it should go forward or backward from the existing node. 

USAGE 

• CALLING SEQUENCE 
CALL SET (IMODE) 

Arguments: 


Parameter 

Name 

In/ Out 

Dimension 

Type 

Description 

IMODE 

In 

1 

I 

Indicator of whether the 
interactive feature is 
needed in SET 


• DATA in/Out 

Labeled COMMON (refer to the labeled COMMON description section): 


Block Name 
Cl 3 


Input 
1 - 501 
502 - 2501 

2502 

2503 


Output 
1 - 501 
502 - 2501 

2502 

2503 
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Block Name 


Output 

2504 


Input 

Cl 3 ^504 

2505 


CIS 

C26 


1-500 1 

501 - 1000 


C-27 1 “ 50 1 

51 

3551 

3651 

4151 

4652 

5152 

5652 


C-28 


7 
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- 500 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

- 50 

• 3550 

• 3650 

■ 4150 

■ 4651 

■ 5151 

- 5651 

- 6252 

1 

2 

3 

4 

5 

6 
7 



Block Name 

Input 

Output 

C28 

8 

8 



9 



10 



11 



16 

C29 

1 - 500 



501 

501 


502 



METHOD 

In each level, subroutine SET determines the subproblem of free variables 
and sends it to subroutine SIMPLX. SET tests the result of the subproblem 
which has been solved by SIMPLX, v.hen directs the program to branch forward 
if there is a possibility of getting a better solution; otherwise, a backward 
step will be taken. 

REFERENCE 

C. E- Lemke, H. M. Salkin, and K. Spielberg, "Set Covering by Single Branch 
Enumeration with Linear Subproblems," Operations Research 19, 998-1022 (1971) 


8/75 
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SUBROUTINE SIMPLX 


IDENTIFICATION 

Name/T i t1 e 
Author/Date 
Machine Identification 
Source Language 


- SIMPLX (Simplex) 

- Han Chang, July 1975 

- UNivAc mo 

- FORTRAN V 


PURPOSE 

Subroutine SIMPLX solves the linear program associated with a subproblem 
defined by subroutine SET- This subroutine provides the LP solution over 
free variables and the information about whether subroutine SET will take a 
forward or backward step. 

USAGE 

• CALLING SEQUENCE 
CALL SIMPLX 

• Data In/Out 

Labeled COMMON (refer to the labeled COMMON description section): 

Block Name Input Output 


Cl 3 

1 - 501 
502 - 2501 

2502 

2503 

2504 

CIS 

1 - 500 

C26 

4 


5 

6 
7 


5/75 
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Input Output 

8 
9 

TO 

1-50. 

51 - 3550 
3551 - 3650 
3651 - 4150 
4151 - 5151 
4652 - 5151 
5152 - 5651 
5652 - 6252 

C28 1 

2 

5 

6 

7 

8 
9 

10 
11 
12 

15 

16 
17 

C29 501 

502 

METHOD 

A revised dual simplex method has been employed in this routine for optimizing 
a linear program. This method was designed to accomplish exactly the same 
function as the original simplex method, but in a way which is more efficient 
for execution on a digital computer. It computes and stores only the information 


1 - 500 
501 


Block Name 


C27 
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that is currently needed, and it carries along the essential data in a compact 
form. In fact, all relevant information at each iteration can be obtained 
immediately after the inverse of basis has been found. The bookkeeping of 
a huge conventional simplex tableau becomes less attractive. 

Another advantage of this method is that the basic inverse is stored in the 
product form rather than a matrix form at each iteration to keep the minimal 
usage of core storage. 

REFERENCE 

G. Hadley, Linear Programming, Addison-Wesley Co,, Inc,, Reading, 
Massachusetts, 1963. 
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SUBROUTINE TABLE 


IDENTIFICATION 

Name/Title . - TABLE CForm Occurrence Table) 

Author/Date - Han Chang» July 1975 

Machine Identification - UNI VAC 1110 

Source Language - FORTRAN V 

PURPOSE 

Subroutine TABLE prints the feasible combination payload occurrence table 

in a particular year by user's request and translates this table into 
input format for subroutine SET, 


USAGE 

9 CALLING SEQUENCE 
CALL TABLE (I YEAR) 


Arguments; 





Parameter 

Name 

In/Out 

Dimension 

Type 

Description 

lYEAR 

In 

1 

I 

The year indicator 
on which data case 
based 


• DATA In/Out 

Labeled COMMON (refer to the labeled COMMON description section); 

Block Name Input Output 

C9 2 

51 

54 


TABLE-1 
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Block Name 


Input 


Output 


C13 1 - 501 

502 - 2501 

2502 

2503 
2505 

Cl 6 1 - 500 

CIS 1 - 500 

501 - 1000 


METHOD 

Subroutine TABLE determines and prints a list for each payload of all fea- 
sible combinations which include that payload by reading the relevant infor- 
mation from a data file {logical unit 2) which has been generated by the MPLS. 


TABLE-2 
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5.4 SAMPLE INPUT/QUTPUT 


This sample input/output 1s to provide the reader with an example of 
executing the SCA interactively. The procedure to sign on the demand terminal 
and execute tLe MPLS is detailed in the SAMPLE User's Guide (.ref, 2 ) and will 

not be repeated here. All the underlined tutorials are the options which the 

user may encounter in the SCA execution. Each of the underlined tutorials 
is accompanied by a section code referring to the location of the option's 
explanation. The alphabet prior to the section code on each underlined 
tutorial will correspond to the short description of the user's response. 

The descriptions are as follows: 

(a) The user wants to execute SCA interactively. 

(b) The user wants unity cost coefficients so he can get a traffic model with 

minimum number of missions. 

(c) The user wants to see the missions with maximum number of payloads. 

(d) The user needs feasible combination number 5 to be included in the 
traffic model . 

(e) The user wants to see the display of Manual Flight/Combination Options. 

(f) The user wants the program to continue to find traffic model. 

(g) After the traffic model is found, the user needs to see the information 
on flight number 7. 

(h) He does not want to see any more, so he enters zero, 

(1) The user wants to have another traffic model, 

(j) He wants flights IQ and 11 to be excluded from the traffic model; the 
missions with maximum number of payloads have been changed. 

(k) The user wants the program to continue to find a traffic model as he 

specified (Notice: flights 10 and 11 did not appear in the traffic 

model . ) 

(l) The user is not interested in seeing the information on traffic model. 

(m) The user does not want another traffic model. 

(n) The user terminates the execution of the SCA. 
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3>RUn HCL:-a.E, ISSOC-iJOU-Ci Fr13-L7^1^7 
date: 09SS75 TIME: 170746 

>5USE TEMPOR. » FD3-L73436*TEMPnP 
RERDY 

>.3MRP TEMPOR. MRP j SAMPLE 
MRP 0026“ 03.^22-17! 09 “<19? > 


STRRT=013533? PRD6 SIZE <I -14225 ^50394 
>Y33*PLIBi. LEVEL 70-1 

BID OF CDLLECTIDM - TIME 3.571 SECOMDS 
: .DXOT SAMPLE 

INPUT TUG CHARACTERISTICS AHD MISSION MODEL DATA: 

<FOR EXAMPLE IPADD SAMPLE. DATR99:* 

>PADD TEMPOR. DATA 15 
SELECT AM DPTIPM: c 3 TO TERMINATE > 

;-i 

SELECT DISPLAY OPTIONS: 7 FOR ALL 5 FDR HONE • 

:.>4 

SELECT AN OPTION: •; 3 TO TERMINATE :< 

>2 

SELECT AN ANALYSIS TYPE: ■•.4 FOP MDNE> 

. k* 

IfiPUT YEAR FDR ANALYSIS: c JO 91 ■ 

•30 

SELECT INTERACTIVE DPTIOHS: fqr r<OME> (a) 3.2. 

T2 ~~ ^ 

• SELECT PERSONAL DATA SRSE TO GENERATE FEASIBLE MISSION: FOR HOME.' 

INPUT MAXIMUM NUMBER OF PR LDATi^ ALLOWED OtfE COMBINATION: 

>3 

SELECT MISSION TYPE: 

1! FOP INPUT CHANGES TO LIST 
O: NO CONSTRAINTS APPLIED 
-l; apply CONSTRAINTS USING LIST 

S-l 

SELECT DISCIPLINE MIX: 

l; FOR INPUT CHANGES tq LIST 
O: NO CONSTRAINTS APPLIED 
“i: APPLY CDNSTRAirJS USU^G LIST 

■•“1 

INPUT 1 TO PRINT MISSION CLASS CODE LISTS QTHERUISE SUP A LlrtE 
.^0 

INPUT 1 TO PRINT PAYLOAD DISCIPLINE MIX LISTS OTHERWISE SKIP A LltC 

> ;i 

MPLS S' ARTED ♦•►♦♦♦^♦♦♦* 

statistical ANALYSIS FDR 1930 
TOTAL NUMBER OF COMBINATIONS GENERATED: 24 

NUMBER OF FEASIBLE COMBINATIONS: 11 

NUMBER OF INFEASIBLE COMBINATIONS: 13 

TOTAL ELAPSED TIME: 47 

<ALL TIMES APE IN MILLISECONDS^ 

AVERAGE TIME PER FEASIBLE COMBINATION: 4 

average time per GENERATED COMBINATION: 1 
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:HDDSE cg;T COEFFICIEHT for ERCH FLISHTt (b) 3.2.2 


1 

1930 

OCCURRENCE TABLE 


PAYLOAD 

MISSIONS 

1,‘ 

1000 

1 9 10 

2> 

1050 

2 11 


2000 

3 3 11 

4.* 

2050 

4 10 

5> 

3001 

5 

o,* 

3002 

6 

7> 

3051 

7 

3> 

3052 

3 


MAX NO. 

SINGLES * 


□CCURREttCE TABLE AMD SCft IMTERFftCE REQUIRED 3>& HILLI3ECDMDS 
CHGDSE criteria POP FLlGHT-^CDMBlHftTION SELECTiaNt <5 FDR HOr<E> fc) 3.2.3 


>t 


YOU WILL GET TUTORIALS IN MODE 1. 

AVAILABLE MISSIONS WITH MR>aMUM NUMBER OF PAYLOADS* 
MISSION PAYLOADS 

3 2 

10 a 

11 2 

1 1 

a i 

CHOOSE MANUAL FLIGHT/CDMEINATIDN OPTION: 

■’S"’"' : 

AVAILABLE MISSIONS WITH MAXIMUM NUMBER OF PAYLOADS* 
MISSION PAYLOADS 

9 2 

10 2 

II . 2 

I 1 

a I 

CHOPSE MANUAL FLISHT/ CaMSINATiaN OPTION; ^ 


fdV 3.2.4 


fe) 3,2.4 


TABLE FOR MANUAL FLISHT/COMBINATION OPTIONS* 

0* LIST OF MANUAL FLIGHT^COMBINITION OPTIONS 
N* ENTER COMBINATION "N" 

-1* NEW SELECTION CRITERIA ARE DESIRED. 

-a* CONTINUE ON TO TERMIMATIDN. 

-3* VIEW ALL COMBINATIONS SPECIFIED SO FAR. 

-4: VIEW INFORMATION ON COMBINATIONS SPECIFIED 
-5* REMOVE LAST SPECIFIED COMBINATION 
CHOOSE MANUAL FLIGHT^^COMSINATIOH OPTION* 

>-Z 


IQ FAR 


ff) 3.2.4 


ALL PAYLOADS CAN BE SCHEDULED IN THE FOLLOWING 
5 G 7 3 

TRAFFIC MODEL GOST IS G 

Total elapsed time in set= 95 

TIME PER MISSION IN MILLISEC.* 15 

HIGHEST LEVEL REACHED » 2 

DO YOU WISH TO SEE INFORMATION ON THESE MISSIONS? 


G MISSIONS. 
10 11 


(q) 3.2.5 


0 * nmw 

-1* PRINT ALL 

-2* PRINT ALL AND SAVE ON SCRATCH FILE 
-3* SAVE ON SCRATCH FILE ONLY 
H * ENTER MISSION “N'* 


5'34 


FLT. MD 


LfiUNCH SITE! ETR 


PRYLQFlDSi LS-Ol 
3091 

SHUTTLE SEQUENCE 
ALTITUDE 
INCLINATION 
TOTAL LENGTH DOWN: 
PAYLOAD MARGINS 64313. 
SHUTTLE DELTRV: 1279. 


LCR A 

3-R 

300. 

23.5 

13. 


TOTAL WEIGHT DOWN: 632.0 


LOAD FACTDRi .01043 


>0 


(h) 3.2.5 

INPUT 



“H"” — 
WHICH 

MISSIONS DO YOU WANT OMITTED ? 



> 1 1 

MO 

>0 

AVAILABLE MISSIONS WITH MAXIMUM NUMBER OF PAYLOABSs 
MISSION PAYLOADS 

3 2 

1 1 

2 I 

3 1 

4 1 

CHOOSE MANUAL FLIGHT/COMBINATION OPTION; 


T ^=2 ^ 

ALL PAYLOADS CAN BE SCHEDULED IN THE FOLLOWING 7 
2 4. 5 o : 

TRAFFIC MODEL COST IS 7 

TOTAL ELAPSED TIME IN SET* 172 

• TIME PER MISSION IN MILLISEC.* 24 

HIGHEST LEVEL REACHED = I 
DO YOU WISH TO SEE INFORMATION ON THESE MISSIONS? 
O; NON E 


Missiors-. 


, {k ) 


U) 3.2.5 


>0 


-l: PRINT ALL 

-2: PRINT ALL AND SAVE ON SCRATCH FILE 
-3; SAVE ON SCRATCH FILE ONLY 
N : ENTER MISSION "14" 

STATISTICS FDR CURRENT FLIGHT SCHEDULE 


AVERAGE NUMBER OF PAYLOADS PEP FLIGHT *1.14 
TOTAL NUMBER OF TUGS REQUIRED * 0 

TOTAL NUMBER OF INITIAL DMS KITS REQUIRED * 0 

TOTAL f4UMBER OF SECOND AND THIRD OMS KITS REQUIRED = 0 

INPUT 1 IF YOU WANT A DIFFERENT SCHEDULE! OTHERWISE SKIP A LINE 


(in) 3.?. 6 


SELECT AH OPTION; < 3 TO TERMINATE..! 


M M-’’ 
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